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Метод ветвей и границ для решения целочисленной 
задачи дробно-линейной оптимизации 


У статті в рамках загальної схеми методу гілок та меж обгрунтовано алгоритм розв'язання задач 
цілочислової оптимізації у випадку дробово-лінійної цільової функції та лінійних додаткових обмежень. 
Ключові слова: дробово-лінійна функція, метод гілок та меж, цілочислове розв'язання. 
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В статье в рамках общей схемьт метода ветвей и границ обоснован алгоритм решения задач целочисленной 
оптимизации в случаєе дробно-линейной целевой функции и линейньїх дополнительньтх ограничений. 
Ключевьт6"е слова: дробно-линейная функция, метод ветвей и границ, 

дцелочисленноге решение. 


Вступ 


Вимога дискретності змінних в явній чи неявній формах зустрічається в багатьох 
практичних задачах оптимізації (1-9). Так, в економіко-математичних моделях, наприк- 
лад, коли необхідно знайти валовий випуск продукції чи кількість поголів'я тварин, 
змінні мають набувати цілих значень. Крім того, побудова саме адекватної реаліям мо- 
делі вимагає поряд з лінійними функціями використовувати дробово-лінійні. 

Аналіз останніх досліджень та публікацій. Для знаходження оптимальних планів 
задач цілочислового програмування застосовують точні та наближені методи |1-16|, серед 
яких досить поширеними є комбінаторні методи. Відомі методи розв'язання задач дро- 
бово-лінійної оптимізації (15), (16), розроблено методи дробово-лінійних задач на 
комбінаторних множинах |4-6|. 
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Виділення невирішених раніше частин загальної проблеми. Однак авторам 
не відомі дослідження задач цілочислової оптимізації з дробово-лінійною цільовою 
функцією. В літературі вони не розглядалися, а тому є доцільним запропонувати по- 
становку таких задач та методи їх розв'язання. 

Мета даної статті - поширити метод послідовного аналізу варіантів в рамках схеми 
методу гілок та меж з використанням на кожному кроці симплекс-методу для оцінки 
розв'язань цілочислових задач оптимізації з дробово-лінійною цільовою функцією, а 
також обгрунтувати алгоритм цього методу, використовуючи ідеї Ленда та Дойга (10). 


Обгрунтування отриманих наукових результатів 


Розглянемо задачу вигляду: знайти впорядковану пару ( У (х є) Бо -) , таку, що 


п 


да сруТО і 
ТО уз тах Гдз тако оо» З загетах /(2) (1) 
2 Фрнува 
« 
заобмежень 5, арх і 57 г б; у бе або. (2) 
і 
г 
жив мезо (3) 
Хі -цілічисла УМіє//,, (4) 


де заіредь Сі д,, аз, Б, - дійсні сталі Міє /,, У є /,, /, позначає 


5977 
множину перших К натуральних чисел П, УЛ 
Для розв'язання (1) - (4) виконаємо наступні перетворення. Спочатку, як і в 
відомих |15) алгоритмах Гоморі чи Дальтона-Ллевеліна, послабимо (1) - (4), відкинувши 
умову (4). Застосуємо до задачі (1) - (3) відображення м, яке задамо співвідношен- 
нями 


1 п 
Хо - чо УззХур У/ є» хек, (5) 


п 
2, ех на, 
41 
де знаменник не обертається в нуль (вважаємо у,» 0 Ух, що задовольняють 
(2), 0)) та перейдемо до задачі з лінійною функцією цілі: знайти 


Е(у )з- тах Е(у)- тах З суне у заго тах Е(у) (6) 
ра 


уєк'Я уєВ'ї! уек'"" 
за умов 
« 
У аг»; - Вухо ерО» мі є» (7) 
уч с) 
уї20 Міє ле 0), М 
5 4уУ, ФУ (9) 


іч 


де у -(Уд» Увез У коза нано, 
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Зауважимо, що /(х )у-Е(у ), де Г(х ) визначається з (1) - (4), а Е(у ) - роз- 
в'язання задачі (6) - (9). 

Задача (6) - (9) є лінійною і може бути розв'язана одним із методів лінійного про- 
грамування, наприклад, модифікованим симплекс-методом. На наступних етапах розв'я- 
зання необхідно використовувати перетворення, обернене до (5), повертаючись до вихід- 
ної задачі та перевіряючи умову (4) для останнього знайденого розв'язання. 

Запишемо алгоритм розв'язання задачі (1) - (4), в основі якого лежить загальна 
схема алгоритму Ленд та Дойг (10). Позначимо к - номер ітерації. Під ітерацією будемо 
розуміти один повний цикл алгоритму. 

1. Відкинути умову (4) і застосувати перетворення (3) до (1) - (3), отримаємо 
(6) -- (9). 

2. Розв'язати лінійну задачу (6) - (9). 

3. Якщо (6) - (9) не має розв'язання, то не має розв'язання (1) - (4), інакше 


нехай ха У - екстремаль задачі (1) - (3). 


53 1) ! ! . 
4. Якщо хо (з пада задовольняє (4), то (109, х) - розв'язання задачі 


(1) - (4), інакше перейти на крок 5. 


б и . ! 5 ' 
5. Визначити найменший ІНДЕКС ) КОМПОНЄНТИ Хі точки х,такої, що Хі - Не 


ціла. 
6. Записати два обмеження, що в області (2), (3) відтинають х: 
хуя (10) 
зак ні, (1) 


де Й - ціла частина й 
7. Застосувати до (10), (11) перетворення (3): 
У; «у | ою 


ул но: 0) 


8. Приєднати до останньої задачі вигляду (6) - (9) обмеження (12) та застосу- 
вати кроки 2 - 4 до розв'язання (6) - (9), (12). Якщо задача (6) - (9), (12) не має 
розв'язання, перейти на крок 9, інакше (Р Су) у) - розв'язання задачі (6) - (9), (12). 

9. Приєднати до останньої задачі вигляду (6)- (9) обмеження (13) та застосувати 
кроки 2 - 4 до розв'язання (6)- (9), (13). Якщо задача (6)- (9), (13) не має роз- 
в'язання, перейти на крок 10, інакше (Р (за), у») - розв'язання задачі (6) -(9), (13). 

10. Якщо жодна із задач вигляду (6) - (9), (12) та (6) - (9), (13) розв'язання не 
має, то задача (1)- (4) теж розв'язання не має у випадку к - 1. Для к»і вибрати для 
подальшого галуження іншу область з вершиною, знайденою на кроці 12 (к -1)4ї 
ітерації, і перейти на крок 4. 

11. Якщо одна із задач вигляду (6)- (9), (12) чи (6) - (9), (13) розв'язання не 


ЧІ 
має, то перейти на крок 4, вважаючи хо у (») , де і - номер точки, що надає 


цільовій функції найбільшого в області Р, значення. 
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12. Якщо обидві задачі вигляду (6) - (9), (12) та (6) - (9), (13) мають розв'язан- 
ня, то для подальшого галуження обрати ту, яка надає цільовій функції більшого зна- 


. -, -- -- 1 5 
чення, 1 перейти на крок 4, вважаючи х - у, (хо) , де і - номер точки у;, що надає 


цільовій функції більше з двох значень БЕ (у, , ) - 1,2. У випадку, коли значення ці- 
льових функцій збігаються, перейти на крок 4 і проаналізувати розв'язання кожної із 
задач. 

Обгрунтуємо вищеописаний алгоритм. Опишемо спосіб галуження, відсікання 
та оцінювання методу гілок та меж (МГМ), виходячи із специфіки розв'язуваної цим 
методом задачі. 

Позначимо Р - допустиму область вихідної задачі, тобто множину точок, що 
задовольняють умовам (1) - (4). Розіб'ємо множину Р на частини, що не мають 


спільних точок, тобто реророрб, рор'ор,зб, де Р - множина 
допустимих розв'язань задачі (1) - (4) при додаванні обмеження (10); р - множина 
допустимих розв'язань задачі (1) - (4) при додаванні обмеження іх і, | «Хі« іх б | |. 


р. - множина допустимих розв'язань задачі (1) - (4) при додаванні обмеження (11). 


Очевидно, що множина розв'язань Др! є порожньою для задачі (1) - (4) і з подаль- 
шого галуження може бути виключена. 

Процес розв'язання можна представити у вигляді дерева, в якому вихідна вер- 
шина відповідає плану х, - оптимальному плану задачі (1) - (3). А кожна зв'язана 
вершина з вихідною вершиною відповідає оптимальному плану такої задачі: знайти 
(1) за обмежень (2), (3) та за додаткового обмеження (10) або (11). 

Кожній з таких вершин (множин Д, допустимих розв'язань відповідної задачі) 


надають оцінку (верхню межу): (0,) т тах / (х) ї Якщо оптимальні плани отрима- 
хер; 


них задач задовольняють умови цілочисельності, то план з максимальною оцінкою 1 
буде оптимальним планом вихідної задачі, інакше необхідно продовжити процес 
розбиття. При цьому кожного разу для подальшого розбиття обирають вершину 
(відповідно множину 0; ) з найбільшою оцінкою, тобто з найбільшим значенням ці- 
льової функції. Для ітерацій к » 1 у випадку, якщо обрана область ДР); не містить точок з 
цілочисловими координатами, використовуємо для подальшого галуження область з 
меншою оцінкою, знайденою на попередній ітерації (крок 10 алгоритму). 


Розбиваючи в процесі розв'язання множину Р) на підмножини р,орер, 
мі 


маємо, що оцінка для будь-якої з них не більша за оцінку для вихідної множини Й, 
тобто для всіх Д; має місце нерівність: 4(0.)«5(Р). 

Правила відсікання (крок 6 алгоритму) враховують той факт, що відсікаються 
тільки ті області Д., які не містять точок, що задовольняють (4). 


Враховуючи спосіб галуження та правила відсікання, має місце наступне твер- 
дження. 

Теорема 1. Алгоритм МГМ, застосований до задачі (1) - (4), знаходить її оп- 
тимальне розв'язання. 


й М ж 
Застосуємо до системи нерівностей, що описують множини 0), В), В», пере- 


творення (5), отримаємо О,, о, 0, . Умова (4) при застосуванні відображення (5) 
набуде вигляду: 
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У; З Худ» (14) 


- 1 
де х, - цілі У є /,, Уд - 


5. ху 4 
іч 
Теорема 2. Обмеження (12), (13) не відсікають допустимих розв'язань задачі 
(6) - (9), (14). 
Доведення. Нехай існує у є 0), така, що у - розв'язання задачі (6) - (9), (14). 
Застосовуючи до у перетворення, обернене до (3), отримаємо точку х, хз у (ус іт 


О скільки задачі (6) - (9), (14) та (1) - (4) еквівалентні, то х - розв'язання задачі (1) - 
(4). Враховуючи, що відображення м задає взаємно однозначну відповідність між 
0 - множиною точок, що задовольняють (1) - (4) та 0 - множиною точок, що 


задовольняють (6)- (9), (14) |5|, маємо 0" - у(0"), р -У 107). а отже, маємо 


хе р'.Однак за побудовою множина Р є порожньою для задачі (1) - (4). Суперечність. 


Таким чином, множина 0! не містить розв'язань задачі (6) - (9), (14). Твердження 


доведено. 
| МОВА п 
Приклад. Оптимізувати функцію ------ 2 ТтаХ за лінійних обмежень 
яр 
ху32Ї2, 2х,- х 59, -х, К4х, 58 таза умов цілочисельності змінних хі, Х». 
Розв'язання. І. Відкинемо умову цілочисельності та застосуємо перетворення (5) 
до «послабленої» задачі, одержимо: 3у,52Уу, З» тах за обмежень -у, -3у, 112, 20, 


2у, 7 У» 79Ус 20, зу, 74Уу, 78у020, у,7У.,21. 


за с А о і -2. 199.13 

2-3.Розв'язавши лінійну задачу, одержимо у 2| - -,--,--|,ЗвідКи ха|--,--|. 
54 18 18 УЗ 
4. Враховуючи, що х, та х, - не цілі, переходимо на крок 5. 
зро дн 
б. Записуємо два обмеження, що відтинають х: х25, х2 6. Застосовуємо до 
обмежень х, 25, х, 2 6 перетворення (5), отримаємо у, - 3), 20, - у, бу, 20. 

7. До лінійної задачі кроку І додаємо обмеження у, - 9у 20 1 розв'язуємо Її. 


ца б оз 5 
Е(у)зозх -Ї 55) 


8. До лінійної задачі кроку І додаємо обмеження - у, "бу, 20 ірозв'язуємо її. 
ЕОд-з хз (6.3). 

9. Враховуючи, що БК (у) -Е (у») з -- переходимо на крок 4 алгоритму. 

Оскільки зей оо) задовольняє умовам цілочисельності, то х, - (6,3) і буде 


зораоб у 8 
оптимальним розв'язанням вихідної задачі. Таким чином, 5 (6,3) - розв'язання. 
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Висновки 


Таким чином, у статті побудовано алгоритм методу гілок та меж (ідейно близький 
до підходу Ленд та Дойг) для розв'язання задач оптимізації дробово-лінійної цільової 
функції з урахуванням умови цілочисельності змінних та лінійних обмежень. 

Перспективи подальших розвідок у даному напрямі. Доцільним видається 
надалі програмно реалізувати даний алгоритм та провести оцінку його складності, а також 
числові експерименти для визначення меж практичного застосування алгоритму. 
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оо АЛ 


О.А. Етеїз, О.А. СПпегпепКо 
Вгапсп апа Воипа Метоа Гог Зоіуїпе те Піеоег РгоБіет 


ої Іпеаг-Етаспопа! Оріїтігатоп 

Тре агисіе соп8ідег5 Бе ргобіега ої орійтігайоп ої ппеаг-Їтасйопа! обуесйує Гапспоп 
оп а роїгупедгоп Ша! 15 Фезсгібед Бу Ше 5узіет ої Ппеаг соп5ігаїпіз апа уліб сопашоп ої Фе 
іпеєег уагіабіез. Асішайу ої Щіз гуре ої ргобіет 15 плойуаїед Бу ргасйїса! пеедз, їп рагії- 
сшаг, ре соп5ігиспоп апа 5їиду ої арргоргіаїе плодєія (о геа! ргосе55е5. 

Тре ашфог5 ргоро5е Фе теоа їог 50Іміпє Ше іпіеєег ргобіет ої Ппеаг-Їтасйопаї 
орйтігацйоп уліб адаїшопаї Ппеаг сопзітатіз. Матеїу, Ше пефШоа ої зедиепна! апаЇузіз ої 
уагіапі5 їп Ше сопіехі ої Ше бгапсп апа бопла пефоай 15 ехіепаеа. Тіз агіїсіе дебсгібеє Фе 
хуаує ої Бгапсріпє, сийпе-оїї апа еуаїшацйоп їп ре Бгапсб апа бойпа пефоа іакіпеє іпіо 
ассоцті Фе з5ресійс5 ої Ше ргобіет Беїпеє 50Іуед. То емаїшцаїе Ше іпіеєег 50Їийоп5 ОЇ ре 
орійтігайоп ргобіет улії а Ппеаг-Їтасйопа! обіесйує Гипсйоп а їгап5ійоп (о а ШПпеаг 
обіесйуєе Гипсйоп аге регіогтеай апа Фе 5ітріех птефоай 15 п5еа оп еуегу 51іер. 

Тре 5оЇийоп аїсогійпт Їог іпіеєег орійпигайоп ул Ше Ппеаг-ігасйопа! обуесйує Гип- 
споп 15 Бий. ТрБе согтесіпе85 ої Ше аїсогіфрта 15 огоппаед Бу Фе Шеогет Шаг Ше ргорозей 
аїєогійрт ої Фе Ббгапсб апа бойла тейоай арріеай (о Фе орійтігайоп ргобіега мії Пе Ппеаг- 
їгасйопа Гипсйоп ул Ппеаг сопяігаїпі5 апа іпіеєег уагіаріе5 Ппда5 Ше орйта! з01ийоп ОЇ 
ргобіет. 

Тре ехатріе ої 50Іміпє Ше ргобіет ої Біз5 гуре бу Ше ргорозеад аїсогіт 15 єїмеп 
їп Фе агіїсіе. 


Стаття надійшла до редакції 25.01.2012. 
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